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Abstract 

An holographic f(R) gravity model of dark energy is proposed. The corre- 
spondence between the f(R) geometrical effective energy density with the holo- 
graphic density, allows the reconstruction of the f(R) gravity in flat FRW back- 
ground in the Einstein frame. The proposed infrared cut-off for the holographic 
energy density depends on two parameters which are fit using the luminosity 
versus redshift data, allowing a suitable reconstruction in two representative 
cases of the EoS parameter: for u> > — 1 and u < — 1. 
PACS: 98.80. -k, 95.36,+x 

1 Introduction 

Many astrophysical data, such as observations of large scale structure |T], searches 
for type la supernovae [2], and measurements of the cosmic microwave background 
anisotropy [3], all indicate that the expansion of the universe is undergoing cosmic 
acceleration at the present time, due to some kind of negative-pressure form of matter 
known as dark energy (for a review see[4j,[5j). Although the cosmological observations 
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suggest that dark energy component is about 2/3 of the total content of the universe, 
the nature of the dark energy as well as its cosmological origin remain unknown at 
present. Among other approaches related with a variety of scalar fields (see [1]), 
a very promising approach to dark energy is related with the modified theories of 
gravity known as f(R) gravity, in which dark energy emerges from the modification 
of geometry [B] , [7J , [E] , [H] , [ID] • Modified gravity gives a natural unification of the early- 
time inflation and late-time acceleration thanks to different role of gravitational terms 
relevant at small and at large curvature and may naturally describe the transition from 
deceleration to acceleration in the cosmological dynamics. Among others, modified 
gravity is one of the areas of the gravitational science in which S. Odintsov has 
contributed to develop thanks to a great number of pioneer investigations. Thus, 
the connection of modified gravity theories with M-theory was suggested in j6j and 
some explicit asymptotic examples were shown. In [7J , [8] , [UJ , [12] , [T3| , a modified 
gravity model with positive and negative powers of curvature and InR terms were 
proposed, to describe different unification scenarios of early time inflation with the 
late-time cosmic acceleration. This unification was also studied with a Gauss-Bonnet 
invariant dependent function as proposed in [H], [15]. The one- loop quantization 
approach to f(R) gravity was investigated in [T6]. The asymptotic behavior of a wide 
variety of f(R) models, for different cosmological epochs in the FRW background has 
been presented in [T7J and a consistency with experimental tests has been studied 
in [T7] , [IB] , [IH] , [20] ■ A comprehensive work on reconstruction of different modified 
f(R) models, has been done trough the works [21] , [22] , [23] , [21] , [22], where it was 
explicitly demonstrated which versions of above theories may be reconstructed from 
the known universe expansion history (see also [10], [2H], [2Z] for f(R) reconstruction). 
Other generalizations like a non-local model of f(R) gravity, non-minimal coupling of 
f(R) gravity with matter, with Maxwell theory and non-minimal Yang-Mills coupling 
were discussed in [281 EH] [30], EI] respectively. The structure of future singularities 
in models of f(R) gravity consistent with local tests, have been worked in [32], [33] . 
and the constraining of some f(R) models according to VIRGO, LIGO and LISA 
experiments have been investigated in [31]. 

Another way to the solution of the dark energy problem within the fundamental theory 
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framework, is related with some facts of the quantum gravity theory, known as the 
holographic principle. Applied to the dark energy issue, this principle establishes an 
infrared cut-off for the so called holographic energy density, related with cosmological 
scales (discussion and references are given in [35j,|36j). An infrared cut-off given by 
the future event horizon was proposed in [35J and a generalized holographic infrared 
cut-off depending on local and non-local quantities was proposed in [37] . Viewing the 
modified f(R) gravity models as an effective description of the underlying theory of 
dark energy, and considering the holographic vacuum energy scenario as pointing in 
the direction of the underlying theory of dark energy, it is interesting to study how 
the f(R) gravity can describe the holographic energy density as an effective theory. 
The reconstruction of f(R) theories under different conditions has been presented in 
[2H E21 E31 E2 ESI EU EJJ and the holographic reconstruction of f(R) gravity have 
been discussed in [27] with the infrared cut-off given by the future event horizon. In 
this contribution, an holographic f(R) reconstruction using the new infrared cut-off 
given in [361 EE], is presented. 



2 The holographic model 

Let us start with the holographic dark energy density given by [361 138] 

Pa = 3Mp (a>H 2 + /3h\ (2.1) 

where H = a/a is the Hubble parameter and at, {3 are constants which must satisfy 
the restrictions imposed by the current observational data. This kind of density 
may appear as the simplest case of more general f(H,H) holographic density in 
the FRW background. Note that this proposal avoids conflict with the causality. 
Here the constants a and (3 will be fixed according to the luminosity versus redshift 
observational data as given in [39], based on [ID] data. As we will focus on late 
time reconstruction, we consider that the holographic dark energy density gives the 
dominant contribution to the Friedmann equation and neglect the contribution from 
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matter and radiation, thus the Friedmann equation becomes simpler 

1 



which gives a power-law solution 



3M p 2 



H 



p A = a H z + (3H 



(3 1 



(2.2) 



a — It 

with constant equation of state parameter given by 

2a - 1 

LdA = -1 + 



3 (3 



In terms of the redshift parameter we can write H in Eq. (2.3) as follows 

H = H (l + z) (a - 1)//3 



(2.3) 



(2.4) 



(2.5) 



where Hq = and to is the present time (z = 0). In the flat FRW background 

the luminosity distance can be written as [1] 

dz 



d I 



1 + z) 



H(z) 



(2.6) 



where we used the light speed c = 1. With H(z) given by (2.5) the luminosity distance 



takes the explicit form from (2.6) 



d L (z) 



1 + z 



P 



1 + z )^~ a+1 )/^ _ i 



(2.7) 



H P-a + 1 

based on this result we can plot the luminosity distance Ld in terms of the redshift, 
and tune the parameters a and (3 in such a way that may fit the data reconstructed 
plot as shown in [S] and [22] • The plot for three representative values of the a and (3 
parameters is given in Fig. 1 
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Figure 1: The luminosity distance \og 10 (HodL(z)) versus redshift. The three curves 
show H dL in logarithmic scale for Q m0 = 0, Q\ = 1, a = 1.1, 1, 0.95 and (3 = 0.5. 
The a = 1 curve corresponds to the de Sitter solution. 



3 Late time reconstruction of f(R) 



Let's start with the action for f(R) gravity with matter 

1 



S 



d X\f—g 



2k 



2 f(R) + £r. 



(3.1; 



where f{R) is a generic function of the Ricci scalar curvature R, and C m is the matter 
fields Lagrangian density. Variation with respect to the metric and setting aside the 
complications with the variations of /(i?)-type actions leads to the field equation 
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f(R)R, u - -f(R)g,u - (V^V, - g^U) f(R) = k 2 T^ 



(3.2) 



where is the covariant derivative, prime denotes the derivative with respect to R 
and T^ is the stress-energy tensor for the standard matter source. Interpreting this 
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field equations in the form of Einstein equations with an effective stress-energy tensor 



composed of curvature terms moved to the right hand side, Eq. (3.2) can be written 
as (hereafter we will take k 2 = 1) 



(m) 



f'(R) ^ 



(3.3) 



with Tp/ given by 



r r(R) = 1 

f'(R) 



(f(R) - Rf'(R)) 9iiv + V„V V /' (R) - g^Uf{R) 



(3.4) 



where T^ v may be defined as an effective stress-energy tensor for some kind of curva- 



ture fluid. This allows to put Eq. ( 3.4 ) in the form of a perfect fluid energy- momentum 
tensor, which will turn out to be useful in our reconstruction scheme. Note also that 
in this case the matter term becomes coupled to the geometry through f'(R). 
Assuming flat FRW background, the modified Friedmann equations take the form 
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3f'(R) [ 



Pm + Rf(R)-f(R) _ 3HRnR) 



(3.5) 



in this way we can define an effective curvature- density 
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PR 



f'(R) l 



Rf'(R) - f(R) 



- 3HRf"(R) 



(3.6) 



The space-space components of (|3.3|) lead to the other Einstein equation 
1 



2H + 3H Z 



— ^ ,p m + Rf"(R) + R 2 f'"(R) + 2HRf"(R) + l - (f(R) - Rf'(R)) 

(3.7) 



which allows the identification of the effective curvature-pressure 
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PR 



f'(R) 



Rf(R) + R 2 f"'(R) + 2HRf"(R) + \ (f(R) - Rf'(R)) 



(3.8) 



Hence, the curvature terms in the r.h.s. of Eqs. (3.6) and (3.8) can be viewed as an 
effective curvature-fluid with equation of state parameter expressed as 

Rf"{R) + R 2 f"'(R) - HRf"(R) 



-1 + 2- 



Rf'(R) - f(R) - QHRf'(R) 



(3.9) 
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defining the function <&(R) = f'{R) [2] and combining Eqs. (3.5) and (3.6) it follows 



that the EoS parameter for the curvature-fluid (3.9) can be written as 



-1 + 



H<5> 



(3.10) 



3# 2 $ - Pm 

where the last term p m is the pressureless dark matter contribution. According to 
the holographic correspondence this EoS parameter should be equal to the one given 



by Eq. (|2.4|). This conduces to the following equation for the function $ 

2a 
~3~ 



$ - - 2 



a 



1 



'Pi 



(3.11) 



P 3/3 

It should be noted here that the cosmological evolution takes place in the Einstein 



frame as is clear from Eq. (3.3), and therefore only in this frame the solution (2.3) 
makes sense as this solution is obtained via solving the standard Einstein equations. 
Changing the variable from cosmic time t to redshift z, the matter density p m becomes 



n „, = 3Hfo m0 (l + z) 3 (3.12) 
here we assumed that p m is conserved separately. After this change of variable, the 



equation (3.11) transforms to 



+ (l + zYH 



dH d<$> 



- 2- 



a 



a 



1 



dz 2 dz dz (3 

The general solution to this equation is given by 



P 



Hfo m0 {l + z) 



(3.13) 



2ptt 



17(0 



'1 



\3-2/p 



9p -6p 2 -2 

+ $ (1 + ^(p-WtP- 1 )^ 8 ?)/ 2 ? 



$+ (1 + z) 



(p-l+V(P-l) 2 +8p)/2p 



(3.14) 



where p = and the constants §q will be determined using local solar system 
consistency. To determine f{R), we first reconstruct the form of f(R) from (3.14) as 
a function of the redshift z. From $(i?) = f'{R) it follows that 



df[R(Z)} 



dR 



(3.15) 



dz dz 

Expressing the scalar curvature R = 6(H + 2H 2 ) in terms of z through the equation 



(3.12) and replacing $ from (3.14), the Eq. (3.15) can be integrated with respect to 
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z, yielding 



/(*) ^C^H^l + zf + Cl{p)Hl<$>+ (1 + 2) ( P +3 +v ^Tj^ )/2p 



+ C 2 (p)H^ (1 + z 



,(p+3- v /(p-l) 2 +8p)/2p 



(3.16) 



where 



8(1 - 2p) 
5p3 _ 9^2 _|_ 2p ' 



cfiv) 



24(2p- 1) 



p(p + 3 ± y/(j> - l) 2 + 8p) 



(3.17) 



writing the variable (1 + z) in terms of R using Eq. (3.12), we can obtain f(R) 
explicitly in terms of R 

3p/2 



f(R) 



p 



6# 2 (2p-l) 



+ 



+ 



P 



m%(2p-i) 
p 

6H 2 (2p-l) 



(p+3+y/(p~l) 2 +8p)/4 



(p+3- v /(p-l) 2 +8p)/4 



Hl^Ct{p)R {p+ ^ [p - 1)2+Sp)/i (3.18) 

^ 2 $0C2~(P)^ (P+3_V/(P ~ 1)2+8P)/4 



The constants of integration <3>q can be determined using a local consistency con- 
ditions with the aim to not affect by measurement errors the reconstructed f{R). 
This conditions are f'{R) z =o — 1 (to recover the actual value of the Newtonian con- 



stant from Eq. (3.3) at z = 0) and f"(R) z=0 = (consistency with the weak field 
approximation), which translates into (see |42j.|26j) 

'd 2 f\ fd 2 R s 



dz 



2 = 



dz J 



2=0 



dz 2 



2=0 



dz 2 



(3.19) 



2=0 



then, using this conditions the constants $q become 



=F 8(2p 2 - hp + 2) =f 4(2p -l)(3-p^ vV + 6p + 1) 
+ p 2 d(p)fi m0 (3 - 5p T Vv^+Gp+l) I 
pC±(p) vV + 6p + lG» + 3 ± vV + 6p + 1) 



(3.20) 



Note that this reconstruction was made in the low redshift region, where the power 
law expansion resulting from the restriction of the holographic model to the dark 



energy dominance (see (2.2)), is an acceptable approximation according to the data, 
as shown in Fig. 1. Note also that our model is approximate as is obtained from 
reconstruction scheme and it is valid within the adopted approximation. It is not so 



easy to recover more realistic non-linear f(R) gravity which gives theory 3.18 as the 
approximation. 



4 Discussion 

In this work we showed a reconstruction of the late time cosmo logical f(R) gravity 
dynamics, based on the holographic energy density cut-off proposed in [36j,|38j. An 
analytical expression for the f(R) model compatible with observational data at low 
redshift has been found. The reconstructed Lagrangian is of a power-law type and 
contains only positive curvature powers for a > l,/3 > (quintessence phase) and 
negative curvature powers for a < 1, (3 > (phantom phase), except in the region 
of parameters —3 — y/8 < j3/(a — 1) < \/8 — 3, which is forbidden in this model 
as the power becomes complex. The a < 1,(3 < can also be considered, giving 



similar to the case a > 1,(3 > results. From equations (2.2) and (2.4) it follows 



that for a = 1 we obtain the de Sitter solution with u\ = — 1. The last two terms 



in (3.18) correspond to the homogeneous solution of (3.11) with Q m0 = 0, which 
describe a purely geometrical-dominance model of dark energy . For the specific case 
of a = 1.1, (3 = 0.5 we expect that the term with lower power R 2 -^/ 4 dominates. 
The matter dominance epoch occurs for p = 2/3, giving the Einstein term plus 
corrections (this corrections may be caused due to the extrapolation to high redshift). 
Is important to note that the holographic dark energy density used here was obtained 
within the Einstein frame and not as the solution to the modified f(R) gravity. Hence 
the reconstructed f(R) theory effectively describes the holographic dark energy in 



Einstein gravity. Nevertheless, Fig. 1 shows that the proposed solution (2.3) is a 
good starting point, at least at the phenomenological level, to reconstruct the f(R) 
model at low redshift independent of the theory. 

In summary, we have studied how the modified f(R) gravity model can be used to 
describe the holographic energy density as an effective theory at low redshift, which 
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conduces to the reconstruction in a direct and unambiguous way. However, this 
reconstruction is settled at the phenomenological level, and the theoretical root of 
the holographic density still to be investigated. 
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